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Abstract 

For any exponential Levy model whose diffusion component is nonzero, we provide an exact series 
representation for the implied volatility of a European call option. Numerical examples are provided. 
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1 Introduction 



Various approaches have been taken to studying the implied volatility surface induced by a given class of 
equity models. Most of these approaches explore asymptotic regimes of strikes and maturities (long ex- 
piries, short expiries, large strikes, small strikes, etc.) or a specific feature of the implied volatility surface, 
such as the at-the-money skew. An exhaustive review of the implied volatility literature would be pro- 
hibitive. But, we mention a few papers that deal with exponentia l Levy processes in particular. The short 



maturity volatility smile is studied in 



Figueroa-Lopez. Forde. and Jacauieij (|201ll) . The model- free results of 



Figueroa-Lopez and Fordd (l2012f). and t h e long matur i ty snr ile in 



tooi ) 



Gao and Led (|2011l ) take 



a particularly simple form for exponential Levy models and, as such, are useful for studying extreme strike 
behavior (large and small) of implied volatility. Fo r a review of re s ults on asymptotics for impli ed volatilities 



in exponential Levy models we refer the reader to 



Tankovl(|201ll) : 



Andersen and LiptonI (|2012l ) 



Our approach to studying implied volatility is quite different from the above-mentioned works. Rather 
than exploit a particular maturity and/or strike regime, we exploit the simple structure of exponential Levy 
models. In doing so, we obtain an exact formula (written as an infinite series) for the implied volatility of a 
given call option. As far as we are aware, this is the first time a formula for the exact implied volatility has 
been given in any framework - exponential Levy or otherwise. We also mention that our formula is extremely 
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simple to derive. While previous authors have used advanced mathematical techniques (e.g., saddle-point 
methods, moment analysis, large-deviation principle, etc.) to derive asymptotic implied volatility results, 
our exact result requires only basic calculus. 

The rest of this paper proceeds as follows. In section [5] we introduce the class of exponential Levy 
models. In section [3l we review how European options may be valued in an exponential Levy setting using 
generalized Fourier transforms. Finally, in section|3]we define implied volatility and - for a given call option 
and exponential Levy model - derive a formula for the corresponding implied volatility. The main result of 
our work is summarized in Theorem [71 Numerical examples are provided at the conclusion of this paper. 



2 Exponential Levy Models 



In this section we review the class of exponen t ial Le vy models. A detailed development can be found in 



Cont and Tankovl (|20041 ) 



0ksendal and SulemI (j2005r ). We assume a frictionless market, no arbitrage and 



take an equivalent martingale measure P chosen by the market on a complete filtered probability space 
{^1,3^, {3^t,t > 0},P). All processes defined below live on this space. Let S represent the price process of a 
risky asset. The main assumptions of this paper is that S can be modeled as an exponential Levy process 



dXt = -fdt + a dWt + / z dNt{dz) 

l\z\<R 



z\>R 



zdNt{dz), 



Xn = X. 



Here, R € [0, oo], the volatility satisfies cr > 0, is a Brownian motion and is a Poisson random measure 
characterized by Levy measure v 



^Nt{dz) = v{dz)dt, 
We require that v satisfy 

min(l, < oo, 



dNt{dz) = dNt{dz) - v{dz) dt. 



I 

Jr 



e^v{dz) < oo. 



(1) 



z\>n. 



The first condition must be satisfied by all Levy measures. The second condition guarantees that ¥. St < cxd 
for all t e M+. Valid choices for R depend on the Levy measure v. We can always choose R — 1. If 
i|2|>i kl^lo''^;) < oo then we may choose R = oo. For simplicity, we assume S pays no dividends and the 
risk-free rate of interest is zero. Thus, S must be a martingale. The martingale condition is satisfied if and 
only if 

1 n 



7 ^ 



iy{dz) (e"^ - 1 - zl{|2|<flj) 



(2) 
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3 European Option Pricing 



We consider a European option expiring at time t > with payoff h{Xt). Using risk- neutral pricing, the 
time-zero value of such an option is the P-expectation of the option payoff 



Lewij (|200ll ): iLiptonI (|2002l ) independently show that u{t,x) can be computed using generalized Fourier 
transforms. We review their method below. For brevity, we do not include any proofs. Let 0(A) denote the 
characteristic exponent of X 

„2 



cpiX) = ijX - yA^ + / i,{dz) (e^^^ - 1 - zAzI{|,|<fl,}) 



</)(A) :=logEe*^^S 
We assume that </> is analytic in an infinite strip A"^ of the complex plane 



A"^ :={A€C:Im(A) € (A^A^)}, 



A!: = inf U < 



I , At = sup |a > 1 : ^ iy{dz)e^' < ooj 



Let h{X) denote the generalized Fourier transform of h{x) 

1 



h{X) := 



dxe-'^'-'h{x). 



We assume h{X) is analytic on an infinite strip of the complex plane of the form A'' := {A e A"^ : Im(A) e 
(A'ijA'j:)}. Let A = A^ + i\ where Xr^K G M and fix the imaginary component: A^ e (A'l,A^j:). Then the 
value of the option u(t,x) is given by 

1 



u{t,x) = / dA^e*'^(^)/I(A)V'A(2^), 



2tx 



4 Implied Volatility 



In this section we fix and a call option payoff h{x) = (e^ — e^)'^ . Note that 



h{X) 



Im(A) < -1. 



/2^(iA + A2) 

We also fix (T > and a Levy measure v — e ji where e > and /i is any Levy measure that satisfies (H]). By 
([5]), the parameter 7 is fixed by a and v. To keep track of their dependence on e we write the characteristic 
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exponent as </'(A) and the option price u{t,x) as 4>^{X) and u^{t,x) respectively. We have 

^^(A) = 0o(A)+e(/<i(A), 
<^o(A)-ia2(-A2-zA), 

<^i(A) = -i\ I fiidz) (e^ - 1 - zl{|^|<j^}) + / fi{dz) (e*^^ - 1 - iXz^^^^Ry) , 

and 

= f (iAe*'^'(^)MA)V'A- (3) 
To ease notation, we have dropped the subscript r from dXr. The fohowing definitions wih be useful: 
Definition 1. The Black- Scholes Price u^^ : R+ R+ is defined as 



(p) := f dAe*^"'(^^'')^(A)VA, </>^^(A;p) = ^p^-X' ^X). 



Definition 2. The Implied Volatility is defined implicitly as the unique number cr^ e such that 
where is given by 

Remark 3. For < i < oo the existence and uniqueness of the implied volatility can be deduced by 
using the general arbitrage bounds for call prices and the monotonicity of u^^ . 

Remark 4. Note that u^^ is an invertible analytic function that satisfies dpU^^{p) > for all p > 0. By 
the Lagrange inversion theorem, the inverse [w^"^]^^ of such a function is also analytic. 

Our goal is to find an explicit formula for the implied volatility . To this end, we note that 

gt0^(A) ^ gt(0o(A)+e0i(A)) ^ gt0o(A) y ]_ (ie^^(A))" . (5) 

Inserting ^ into ^ we obtain the following series representation Q for u'^ 

?x^ = Ve"u„, 7.„ = - / dAe*^°(^)(0i(A))"/I(A)V'A. (6) 

Note in particular that uo = u^^{a). 

From (l6|), it is clear that is an ana lytic function of £. I t is a useful fact that the composition of 



two analytic functions is also analytic (see 



Brown and Churchilll (|l996l ). section 24, p. 74). Thus, in light 



^By Fubini's Theorem, exchanging the order of summation and integration is allowed since /jj le**^ ^■^"> h{\)4>\ \ < oo. 



of Remark 21 we deduce that = [u^^] ^{u^) is an analytic function and therefore has a power series 
expansion in e. We write this expansion as foUows 



a" = (70 + (5% 
Taylor expanding u^^ about the point ctq we have 

OO ^ 



(7) 



fc=i 



n=0 
„.BS 



OO / OO \ " 

OO 
OO 

^"M + E^' 



E E n-. 



fc=l 



(^o)+E^' 



OO _j / n \ 

-^^^+E;^ E n-.J'" 

n=2 \jl+■■■+],^=ki=l I 



(8) 



Now, we insert expansions (|6]) and ([8]) into (|4]) and collect terms of like order in e 
0(1) : "0 = u'^Hgo), 



Uk 



«=2 \ji+---+j„=fci=l / 



fc > 1. 



Solving the above equations for {cfe}^Q we find 



0(1) : <Jo = <J, 



0(e'=) : dfc 



fc > 1. (9) 



Remark 5. The right hand side of ([9]) involves only aj for j < fc — 1. Thus, the {crk}'^i can be found 
recursively. 

Remark 6. Note that d^u^^{a) is easily computed using 



a>^^(a) = fdX (9^e*^°(^)) /I(A)^A. 
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Explicitly, up to 0(£^) we have 
0(e): 

0{e') : 
0{e') ■■ 



ai = 



0-2 



(^3 = 



(74 



Ul 

^3 - (q-2gi'9g + jj(jfdl)uo 

daUo 

M4 - ((73(71 + ^(7pg + l(J2(^ldl + ^(7fl9^)Mo 



We summarize our main result in the following theorem: 

Theorem 7. The implied volatility (7^ defined in (j4]) is given explicitly by ^ where ctq = (7 and {(7fe}^i 
are given by ([9|). 

Remark 8. We emphasize: we have made no assumption about the size of e. Theorem [7] is valid for any 
e > 0. In particular, one can always choose e — 1. 

Remark 9. Everything we have done so far is exact . The accuracy of the implied volatility expansion ([7]) 
is limited only by the number of terms one wishes to compute. 

Define the 0(e") approximation of the implied volatility 



E 

/£ = 



At the end of this document, we provide numerical examples illustrating convergence of cr^"-* to for three 
well-known exponential Levy models: 

• the Jump-diffusion model of [ 



MertonI (|1976l) : figure [H 



• the Variance Gamma model of 

• the CGMY model of 



Madan. Carr. and Chand |l998l ): figure H 



Carr. Geman. Madan. and Yoij (|2002i ) : figure |31 



We plot implied volatility as a function of the log-moneyness to maturity ratio, LMMR := [k — x)/t. In all 
three models, we see excellent convergence of cr^"^ to (7^. Convergence is fastest for values of k near x and 
slows as k moves away from x. 

Remark 10. Although our focus has been on exponential Levy models, the exact implied volatility expansion 
outlined above will work for any model whose European call price can be expanded analytically in e as 

oo 
k=l 



where e is some model-specific parameter. See, for example, 



Lorid (12012X1 1 
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